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We describe a calculation of heavy-light decay constants including virtual quark loop effects. We have generated 
dynamical gauge configurations at three (3 values using two flavors of Kogut-Susskind quarks with a range of 
masses. These are analyzed with a Wilson valence quark action. Preliminary results based on a "fat-link" clover 
valence quark action are also reported. Results from the two methods differ by 30 to 50 MeV, which is presumably 
due to significant — but as yet unobserved — lattice spacing dependence in one or both of the approaches. 



Decay constants for the B and Bg mesons 
are crucial for the accurate determination of the 
CKM mixing matrix. Reference 0| describes 
our evaluation of these decay constants in the 
quenched approximation; the results are consis- 
tent with those from several other groups ||] . The 
effects of quenching in ^ were estimated by com- 
paring with results including dynamical quark 
effects at fixed lattice spacing. We now have 
enough results with Np = 2 dynamical quarks 
to start to study the continuum limit in the dy- 
namical theory. This is the crucial step to go from 
quenched answers with estimates of quenching ef- 
fects to true dynamical answers. 

Dynamical gauge configurations have been gen- 
erated with two flavors of staggered quarks at 
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(3 — 5.445, 5.5, and 5.6, with a range of dynami- 
cal masses. (See Table 1^.) We have analyzed each 
set with Wilson valence quarks (both heavy and 
light) as well as static heavy quarks, as in ||l[. In 
addition, we have begun to use heavy and light 
"fat-link" |^,^ clover valence quarks. 

In the fat-link clover case, we implement the 
full Fermilab program Q through 0{a) and 
through 0(1/M), including the 3-dimensional ro- 
tations ("di" terms). The shift to the kinetic 
mass is done as in the Wilson case [0, except 
that tadpole improvement is not needed. 

In general, we treat the dynamical quark con- 
figurations as fixed backgrounds and perform chi- 
ral extrapolations in the valence quark mass only; 
i.e., we do "partial quenching." However, we 

also try extrapolating with mvalencc = "Idynamical 

("full unquenching" ) . The difference is treated as 
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Table 1 

Lattice parameters. All sets use Np — 2 dynam- 
ical staggered quarks and are analyzed with Wil- 
son valence quarks. To date, 98 configurations of 
set R have been analyzed with fat-link clover va- 
lence quarks. Set G was generated by HEMCGC. 



name 




aniq 


size 


# configs. 


L 


5.445 


0.025 


16^ X 48 


100 


N 


5.5 


0.1 


243 X 64 


100 





5.5 


0.05 


243 X 64 


100 


M 


5.5 


0.025 


20^ X 64 


199 


P 


5.5 


0.0125 


20^ X 64 


199 


U 


5.6 


0.08 


243 X 64 


201 


T 


5.6 


0.04 


243 X 64 


202 


S 


5.6 


0.02 


243 X 64 


201 


G 


5.6 


0.01 


16^ X 32 


200 


R 


5.6 


0.01 


243 X 64 


200 



a systematic error, although in most cases it is 
smaller than the statistical errors. 

Most other systematic errors (excited states, 
chiral extrapolation, fitting errors in 1/M, pertur- 
bation theory [in Wilson case] , difference between 
m2 amd [in Wilson case]) are estimated the 
same way as in the quenched approximation jl]. 

Finite volume errors are estimated by compar- 
ing results of sets G and R. Since set G has a 
smaller physical volume than all other runs, this 
is an overestimate. 

The hardest errors to control with our data are 
discretization errors and the effects of omitting 
the dynamical strange quark. We discuss them 
below. 

Details about the fat-link clover approach can 
be found in Refs. ^J^. Throughout the current 
work, we use iV = 10 smearing steps and smearing 
parameter c = 0.45 (c/6 is the coefficient of the 
staple sum). This amount of fattening completely 
suppresses exceptional configurations in the range 
of masses we are studying [^ . With the standard 
("thin link") nonperturbative clover action, we 
found exceptional configurations to be a very se- 
rious problem on our perforce somewhat coarse 
dynamical lattices. 

The clover coefficient has been chosen equal to 
the tree- level value, Csw = 1- The fact that 



fattening suppresses perturbative corrections 
leads us to expect that this value should be very 
close to the all-orders (in g) value for our fat links. 
We plan a nonperturbative evaluation to check 
this. 

Bernard and DeGrand [Q] have computed fat- 
link clover Z factors in perturbation theory. For 
light- light (11) and static-light (si) Za, they find: 

9' 
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(-0.241) 
2 (3 log (aMs) + 0.393) 
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q* = 0.71/a. For Z%, q* has not yet 
been calculated; we use the light-light q* . The 
mass-dependent heavy-light Z^ has also not yet 
been computed. We expect that for moderately 
large masses, the difference between Z^ and Z^ 
will be small: such finite numbers are strongly 
suppressed by fattening. We currently use Z^ for 
heavy-lights. 

At present, we have analyzed only a subset of 
one lattice set (R) with fat-link clover valence 
quarks. With only two light quark masses cur- 
rently available, we choose to focus here on Jb^ ■ 
Figure |l| shows /b, as a function of a in both 
the Wilson and fat-link clover cases. The Wilson 
valence points are consistent with constant be- 
havior in a; allowing a linear term in the fit makes 
almost no difference in the extrapolated value at 
a = 0. However the extrapolated values are in- 
consistent with the fat-link clover result. Possible 
explanations for this discrepancy are: 

(1) The apparently constant behavior of the 
Wilson results is misleading. Indeed one expects 
the Wilson results on dynamical configurations 
to decrease as a — > with roughly the same slope 
as in the quenched Wilson case. (The quench- 
ing effect on this slope should be roughly like 
the quenching effect on physical quantities, i.e., 
~ 5-30%.) In this scenario, the reason that the 
Wilson results look constant is that the effects of 
dynamical quarks (which should raise decay con- 
stants by deepening the potential well at r = 0) 
are turning on as the lattice spacing becomes fine 
enough to see the small r behavior. For smaller 
a, they would begin to fall. If we assume that 
the a = limit of the Wilson data is equal to the 
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Figure 1. /b, vs. lattice spacing. 

fat-link clover result, and that the linear slope is 
the same as in the quenched case, we can make 
a quadratic fit to the Wilson data with a con- 
fidence level of 0.23 and a reasonable quadratic 
term of scale (390 MeV)^. This does not validate 
the scenario, of course, but only shows that it is 
a consistent possibility. 

(2) Too much fattening has done violence to 
the physics governing Jb^ ■ This could be the case 
if, e.g., the smearing softens the Coulomb poten- 
tial at the origin enough to reduce significantly 
the decay constants. This would not mean that 
fattening is "wrong," but that this much fatten- 
ing introduces significant lattice spacing depen- 
dence. This dependence presumably would occur 
at 0{a^) or higher, since we have argued that 
Csw = 1 is close to the nonperturbative value 
needed for 0{a) improvement. 

(3) Perturbation theory, used to find the renor- 
malization constants in the fat-link case, has bro- 
ken down. This may be the case because the fat- 
tening has so reduced the large q behavior of the 
integrands that the integrals are IR dominated, 
and the resulting effective coupling constant is 
too large. The small values obtained for q* are in- 
dicative of this potential problem. In retrospect, 
less fattening would have been preferable 0. 

The full explanation is probably some combina- 
tion of these three scenarios. Scenario (1) makes 
it clear that, while extrapolating the Wilson va- 
lence results with a constant may produce signif- 
icant systematic errors, it should give an upper 



bound to the correct result. If scenario (2) were 
the only problem with the fat-clover data, then 
the fat-clover result would be a lower bound to the 
correct result. Scenario (3) complicates the situa- 
tion. However it is unlikely that the correct result 
is much below the fat-clover result, because it is 
unlikely that the Wilson data would have a slope 
much steeper than in the quenched case. Thus we 
average the constant-extrapolated Wilson and the 
fat-clover results, and use the spread to estimate 
the discretization error. Clearly, this analysis is 
preliminary; much more study is needed. 

With Np = 2, we are missing the effect of a 
dynamical strange quark. To estimate this effect, 
we assume that each dynamical quark, indepen- 
dent of its mass, has the same effect on the decay 
constants. This assumption is supported by the 
Wilson valence data. The values of /_B,in Fig. |^, 
e.g., do not depend strongly on the dynamical 
quark mass (which varies from ^ms/2 to ^4ms). 
We thus estimate the effect of the missing strange 
quark by taking 50% of the difference between the 
Np ~ 2 results and our older quenched results [Q . 

With the above caveats, our preliminary results 
are (in MeV for the decay constants): 
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/b = 194(3)(22)(+^' 
/z5 = 211(2)(27)( 



+ 10n 
) 



=1.12(1)(5)(+ 



:235(2)(if,)(ii3) 



^=l.ll(0)(t*)(ti) . 
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The errors are statistical, systematic (within 
the Np = 2 "world"), and systematic (due to the 
missing dynamical strange quark), respectively. 
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